CalcOpt();

dif0 = diffParamMatrix(0, 10000);
%dif15_10 = diffParamMatrix(1.5e+6, 10);

%dev15 = deviationMatrix(1.5e+6, 365);
%dev75 = deviationMatrix(7.5e+6, 365);
%dev150 = deviationMatrix(15e+6, 365);

dev15log10 = deviationMatrixlog10(1.5e+6, 365);
dev75log10 = deviationMatrixlog10(7.5e+6, 365);
dev150log10 = deviationMatrixlog10(15e+6, 365);

cols0 = ledta(dif0,0.01);
cols15 = ledta(dev15log10,0.013);
cols75 = ledta(dev75log10,0.1);
cols150 = ledta(dev150log10,0.1);

colsToChange = intersect(intersect(cols0, cols15), cols75);

save('difdev.mat', 'dif0', 'dev15log10', 'dev75log10', 'dev150log10');

Prepare;